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The Rabi model describes a two-level system coupled 
with a single bosonic mode [1]. Such a simplest interact- 
ing quantum model has possessed wide applications in 
many fields of physics [2-4]. Therefore, in order to un- 
derstand thoroughly the physical properties of these sys- 
tems, its exact solution, which has been a long-standing 
problem, is of special significance. However, Braak p- 
resented an analytical solution of the Rabi model sev- 
en years ago [5]. The energy spectrum consists of two 
parts, i. e. the regular and the exceptional spectrum. 
Such a spectrum structure is considerably strange. In 
this Comment, I point out that the Braak’s analytical 
solution doesn’t exhibit the true energy spectrum of the 
Rabi model due to the derivation error in solving the 
time-independent Schrodinger equation. 

Braak started from the Rabi model 


Hp =wata+goz(at +a) + Ace (1) 


in Ref. [5]. After taking the transformations a > i and 
at — z, then the Hamiltonian (1) becomes 


_ ( wz; +g(z + 02) A 
Asp = ( A wz0, —g(z+0z) J ` (2) 


Suppose that (ġı $2)” is the two-component wave func- 
tion of Hs. Then one has a coupled system of the first- 
order differential equations 


(z+ 9) gebr(z) + (gz — E)ġi (z) + Ada(z) =0, (3) 
(z — 9) $ 2(z) — (gz + E)bo(z) + Agi(z) = 0, (4) 


where w = 1 and E is the corresponding eigenvalue. 
Braak found that Eqs. (3) and (4) have the following 
solution 


pile) = eF ES o Ky (x) ASH", (5) 
plz) = e EX y Kn(x)(z + 9)”, (6) 


where x = E+g?, E can take an arbitrary value, and the 
constants Kn(x) satisfy the recursive relation (4) in Ref. 
[5]. Obviously, ¢1(z) and $2(z) are divergent at z + —oo. 
Therefore, this two-component solution (¢1 ¢2)7 of Hsp 
is trivial and non-physical due to the divergence of the 
wave function and the undetermined eigenvalue. 

In order to fix the eigenvalue E, Braak employed the 
unitary transformation 


ARS) o 


where the operator T satisfies T(f)(z) = f(—z). It is 
easy to get 


H, 0 
+ P + 
Ut Ha = ( i cae (8) 


Here, H+ = wz0, + g(z + -) + AT. Obviously, H_ can 
be obtained from H, by letting A be —A. 
The time-independent Schrodinger equation for H+ 
with positive parity reads 
d 


aile) +E + We) = Eple) - Apie), 9) 


which becomes 


2ft) — 9 + 2)¥(-2) = BY(—z) — Ay) (10) 


after manipulating T on two sides of Eq. (9). Here w = 1 
and E is the eigenvalue of H} (or Hs»). It is obvious 
that ¢(z) and w(—z) = T(w)(z) in Eq. (9) or (10) are 
correlated due to the reflection operator T. 

With the notation (z) = ¢1(z) and w(—z) = ¢2(z), 
Eqs. (9) and (10) lead to Eqs. (3) and (4), respectively. 
Such a notation is the solution of the coupled equations 
(3) and (4) rather than the single equation (9) with the 
presence of T. I note that if and only if 


G+(x; 2) = $2(z) — T1 (2) = ¢2(z) — g1(-z) =0 (11) 


for any z, this notation {¢(z), ~(—z)} is the solution of 
Eq. (9). 

Obviously, Braak treated (z) and w(—z) as indepen- 
dent wave functions and neglected the condition (11). By 
requiring the wave function {w(z),~(—z)} to be contin- 
uous at z = 0, i.e. G(x; z = 0) = 0, Braak obtained the 
eigenvalues FẸ (see (3) and Fig. 1 in Ref. [5]). Howev- 
er, the constraint (11) does not hold for nonzero z under 
these eigenvalues E [see the expressions (5) and (6)]. So 
such a wave function {w(z),~(—z)} with a cusp at z = 0 
and the corresponding eigenvalue F are not these of H}. 
Similarly, the solution for H_ with negative parity can 
be obtained by replacing A with —A. Therefore, it is out 
of question that the energy spectrum shown in Figs. 2 
and 3 in Ref. [5] is not that of the Rabi model (1) [6]. 

Braak also applied the similar technique to the gener- 
alized Hamiltonian by adding ec, to (1) (i.e. (7) in Ref. 
[5]). However, the same derivation error occurs. The 
energy spectrum depicted in Fig. 4 in Ref. [5] is also 
incorrect [6]. 
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Finally, I would like to mention that I have exactly 
solved the Rabi model by using the unitary transforma- 
tion technique in the occupation number representation 
[7]. The complete energy spectrum is comprised of two 
double-fold degenerate sub-energy spectrum I and II. 
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